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Abstract 

The Casimir force between two parallel plates separated by anisotropic media is investigated. We 
theoretically calculate the Casimir force between two parallel plates when the interspace between 
the plates is filled with anisotropic media. Our result shows that the anisotropy of the material 
between the plates can significantly affect the Casimir force, especially the direction of the force. 
If ignoring the anisotropy of the in-between material makes the force to be repulsive (attractive), 
by contrast, taking the anisotropy into account may produce an extra attractive (repulsive) force. 
The physical explanation for this phenomenon is also discussed. 
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I. INTRODUCTION 


The Casimir effect originating form the quantum fluctuations of the electromagnetic field 
is one of the most remarkable macroscopic effects of quantum physics [l Q . Great effort 
has been put into both theoretical and experimental studies on the Casimir effect [ 4 , |j], as 
it plays an important role in various fields of physics [l]. 

As we know, Lifshitz’s original theory for the Casimir force is only applicable to isotropic 
dielectrics [ 2 , 3]. However, anisotropy can bring us new features about the Casimir effect, 
like the Casimir torque. After the early works on the van der Waals interaction between 


anisotropic bodies in 1970s 6, 3], Munday et al numerically calculated the Casimir torque 
between parallel birefringent plates immersed in liquid, and proposed an experiment to 
observe this torque [8f. As supplemental works, Shao et al gave the analytical expressions of 
the Casimir torque and repulsive Casimir force between two birefringent plates with constant 
permittivity and permeability [ 9 ]. In one of our previous works, we calculated the Casimir 
torque between two parallel anisotropic plates with nontrivial permeabilities, and discussed 
the impact of magnetic properties of the plates on the torque jioj]. On the other hand, 
anisotropy can also affect the Casimir force directly. Romanowsky show that the orientation 
of the optical axis could have significant influence on Casimir force between highly anisotropic 
11]. I11 one of our previous works, we calculated the Casimir force between anisotropic 


plates 

metamaterial plates [12] . The results show that the direction of Casimir force could change 
with both the separation and the anisotropy. Ran Zeng and his collaborators investigated 
the Casimir force between anisotropic single-negative metamaterial slabs and show that the 
electromagnetic responses of the metamaterial parallel and perpendicular to the optical axis 
affected the Casimir force differently [3]- Recent years, the research works have extended 
to the interaction between anisotropic particles and a surface 14, Il5 ]. 

Most of the previous works 7pl|| were based on the assumption that the region between 
the two boundaries was vacuum or filled with isotropic media. However, in fact, this re¬ 
gion can also be anisotropic. Some liquid, such as nitrobenzene and liquid crystal, can be 
anisotropic under special circumstances. One might naturally ask: what new phenomena 
can be seen if the region between the two slabs is filled with anisotropic media? Parsegian 
first calculated the van der Waals energy between two anisotropic bodies acting across a 
planar slab filled with a third anisotropic material in the non-retarded case jf|. Kornilovitch 
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investigated the Van der Waals interactions between flat surfaces in uniaxial anisotropic 
media in the non-retarded limit and discussed the effect of nonzero tilt between the optical 


axis and the surface normal on the interaction 1161. On the other hand, repulsive Casimir 

in 

force is always of great interests to the researchers 2-5, 8 , 9, [121, |13f. According to Lifshitz’s 
theory, when the space between the two slabs is filled with isotropic media, the Casimir force 
can be repulsive if the permittivities of the plates (ei and 62 ) and the interspatial media 
(63 ) satisfy the relation ei < 63 < 62 (or ei > 63 > 62 ) 2, 3]. However, if the interspatial 
media is bircfringent, £3 will have different values in different directions. Naturally, another 
question might be raised: when will the Casimir force be repulsive (or attractive) in this 
situation? Unfortunately, no clear clarification to this has been reported. Without doubt, 
the two problems above are of great importance. Therefore, it is necessary to study the 
possible new features of the Casimir effect between parallel plates separated by anisotropic 
media. In this work we use the quantized surface mode technique 9|, lltj, LUl, [lT] to calculate 
the Casimir force between two isotropic plates when the interspace between them is filled 
with anisotropic media. Our major concern is focused on the impact of the anisotropy on 
the Casimir force, especially the direction of the force. The result shows that the direction 
of the Casimir force can change with the anisotropy of the interspatial media. The physi¬ 
cal understanding of this phenomenon is also investigated. The detailed discussion will be 
presented in the following sections. 


II. THE CASIMIR FORCE BETWEEN PARALLEL PLATES SEPARATED BY 
ANISOTROPIC MEDIA 

The system considered is shown in Fig. 1. Two isotropic plates (with the diameter D and 
the thickness d) made of different materials are kept parallel to each other and separated 
with a distance a. The region between the plates is filled with uniaxial media. The plates 
and the media between them are considered to be nonmagnetic. The x-y plane is chosen 
to be parallel to the surfaces of the plates. It should be noted that if we only consider the 
anisotropy of the media between the plates, but consider the plates to be isotropic, there 
will be no Casimir torque. I 11 this work, the optical axis of the anisotropic media is chosen 
to be in 2 direction, which is often referred as the out-of-plane case. 

Assuming D and d are much greater than a, it is reasonable to disregard the edge effect 
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FIG. 1. Two isotropic plates are separated by anisotropic media with a distance a. The 
surfaces of the plates are parallel to the x-y plane. The optic axis of the interspatial 
anisotropic media is in z direction. 


and finite thickness effect. And the space can be approximately considered to be divided 
into three regions with corresponding permittivities, as shown in Fig.2. The relative per¬ 
mittivities can be expressed as diagonal matrixes, respectively, as following. 


G = 


^ 61 0 0 ^ 

0 Cl 0 
0 0 6 ! 

/ e 2 0 0 ^ 


£2 = 


V 


0 e 2 0 
0 0 e 2 


£3 = 


( e Zx 0 0 ^ 

0 633, 0 

y 0 0 £32 y 


( 1 ) 

( 2 ) 

(3) 


where the subscripts z and x indicate the components parallel and perpendicular to the 
optical axis, respectively. 


According to the quantized surface mode technique mm fl ], we only need to consider 
the zero point energy associated with the surface modes q which are exponentially decaying 
when z>a and z<0. Because each region in Fig.2 is considered as homogenous, the electric 
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FIG. 2. The schematic configuration. The space can be approximately considered to be 
divided into three regions. The optical property of each region is described by 
corresponding relative permittivity. 


and magnetic fields of the surface mode q can be expressed as 


, id m y, is 


E q = iN[a q e q ( k) - aje*(k)], 

(4) 

H (/ = fV[a g hq(k) + ajh*(k)], 

(5) 


where N is the normalization factor, and the parameters a q and a) q are the usual creation 
and annihilation operators, respectively. The parameters e q and h g are the electric and 
magnetic field polarization vectors, k = (. k x , k y , k z ) is the wave vector. We can choose the 
wave vector in x-y plane to be parallel to the x direction, and then the wave vector can be 
written as k = (k x , 0, k z ) = K 0 (a, 0, 7 ), with K 0 = cuc~ x . 

We introduce M\ = €i/e 3x , M 2 = e 2 /e 3x and M 3 = e 3z /e 3x to describe the relative values 
of ei, e 2 and e 3z to e 3x . M 3 can describe the anisotropy of the material filled between the 
plates. M 3 = 1 refers to isotropic. The Casimir energy per unit area at zero temperature 
can be expressed as (£ = — icu) 

fc /*oo poo 

= ^ 2^2 J P d 'P J e 3,xC 2 d^[lnGi(iO + lnG 2 (iO] (6) 

where c is the speed of light in vacuum, and h is Plank constant divided by 27 t. The variable 
p is introduced as p 2 = 1 — a 2 /e 3x . The detailed representations for the functions G*i(i£) 
and Gr 2 (i£) in Eq.(6) are expressed as (The detailed derivation of Eqs. (ED-dS]) is presented in 
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Appendix A.) 


nres_i (s 1 -p)(s 2 -p) 
l 0 (»1+P)(S2+P) 

G2(i0=1 -^:^^fa : ^; exp 


2 pat; 
c 

2 Pat; 


V^3x , 


c ^) ’ 


(7) 

( 8 ) 


(s 1 + M 1 P)(s 2 + M 2 P) 

where s 1)2 = \J M^o. — 1 + p 2 and P = y/ (M 3 — 1 + p 2 )/M 3 . The Casimir force on the plates 


per unit area is 


F = 


P- 


dE(a) 

da 

1 -G 1 

~gT~ 


h 


+ p 


27T 2 C 3 .7! 

1-G 2 “ 

Go 


pdp / «lt 2 ? 3d € 


(9) 


As G 2 is a function of M 3 , the Casimir force F will also depend on M 3 which refers to the 
anisotropy of the material between the plates. And it can be found that G\ does not depend 
on M 3 , which means the anisotropy affects the force mainly through the second term in the 
brackets in Eq. (J9J) . 

For the case that the material between the plates is isotropic, 63 ^ = 63 . = 63 (M 3 = 1 and 
P = p), Eq. (ED becomes 


^isotropic _ 


h 


2i r 2 c 3 
S\ + p s 2 + p 


p 2 dp 


3/2 A 


-exp 


/ 2 paj 


+ 


s 1 -ps 2 -p 

si + Mip s 2 + M 2 p 

-exp — 

s 1 - M\p s 2 — M 2 p V c 


vi)- 

{ 2 pa£ 


3 m 

-1 


Ve 3 



( 10 ) 


which recovers Lifshitz’s result about force on two bodies separated by a gap filled with 
a third isotropic media (equation 4.14 in Ref. Q). 

Let’s turn to the limiting case that the separation a is larger than the characteristic 
absorption wavelength of the material. In this case, we can replace e 1 , e 2 , e$ x and e^ z by 
their values at £ = 0, i.e. the static dielectric constants [3]. The approximate Casimir force 
can be expressed as following 


P ~ 


3 he 

167r 2 a 4 y/ete 


( 11 ) 
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with 


*{M u M 2 ,M 3 ) 



f s w -p \ 

\Sw+p) 


f S 20 — p\ 1 / S\o — MiP\ 

\s 20 +p) pi\s 10 + M 1 p) 


f S 20 — M 2 P \ 

VS 20 + m 2 p ) 


( 12 ) 

where si 0 and s 20 are the values of sx and s 2 at £ = 0. (The detail of the approximate 
calculation of the integration in the Casimir force function is presented in Appendix B) 


III. THE IMPACT OF THE ANISOTROPY OF THE MEDIA BETWEEN THE 
PLATES ON THE DIRECTION OF THE CASIMIR FORCE 

The positive value of F (or T) in Eq. (TITj) corresponds to the attractive force, while the 
negative value corresponds to the repulsive force. As we are interested in the impact of the 
anisotropy of the interspatial media on the direction of the Casimir force, we can just discuss 
the sign of the function T. From Eq. flTTj) and Eq. (fT 2 j) . it is clear that if the permittivities 
satisfy the relations e\ < e 3a; < e 2 and e± < e 32 < e 2 (or e\ > e 3a; > e 2 and e± > e 3z > e 2 ) at 
the same time, T will be minus and the force will be repulsive. This is because although 63 
has different values in different directions, all the possible values of e 3 are still in the range 
form e 2 to . This can make the force repulsive, which is similar to the Lifshitz’s result [3]. 
However, what we mainly concern is not this case but the case when e 3x or e 32 is out of the 
range (e 2 ,ei) (Here, (e 2 ,ei) means the range form e 2 to ei). 

First, we come to the case that only e 3a; is in the range (e 2 ,ei). We let e\ > e$ x > e 2 
(Mi > 1 and M 2 < 1). Fig.3 shows the value of T for different M 3 . The positive value of 
T corresponds to the attractive force, while the negative value corresponds to the repulsive 
force. Two vertical dash-dot lines, corresponding to M 3 = M 2 and M 3 = M \, divide Fig.3 
into 3 regions. I 11 the middle region, where M 2 < M 3 < Mi (e 2 < e 3z < 6i), the force is 
always repulsive, which has been discussed previously. The other two regions (left and right) 
are the regions that we are interested in. In the left region, where M 3 < M 2 (e 3z < e 2 ), 
the force is not always in the same direction and it can be either repulsive or attractive 
depending on the value of M 3 which refers to the anisotropy. In the region on the right, 
where M 3 > Mi (e 3z > ei), the result is similar and the force can also be either repulsive or 
attractive depending on the value of M 3 . 
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FIG. 3. 'IF 'F 1 and T 2 vs. M 3 . Mi = 1.5 and M 2 = 0.8. The two vertical dash-dot lines 
divide the range of M 3 into 3 parts corresponding to e 32 < e 2 , e 2 < e 3 ~ < ei, and e 32 > ei 
form left to right, respectively. The positive value corresponds to the attractive force, 
while the negative value corresponds to the repulsive force. 


In order to see this more clearly, we can resolve T in Eq. (fl2]l into two parts as flowing 


T(M!,M 2 ,M 3 ) = + t 2 


with 


Ti(Mi,M 2 )= / pdp 


1 sip-p s 20 - p 
_P 3 sio + PS 20 +P. 


T 2 (Mi,M 2 ,M 3 ) = / pdp 


1 sio — M W P s 20 — M 20 P 
P 3 Sio + MiqP s 20 + M 20 P 


(13) 

(14) 

(15) 


It is clear that dF is independent of e 32 and M 3 , which means anisotropy of the in-between 
media does not affect \Ef. If > e 3x > e 2 (Mi > 1 and M 2 < 1), as we assumed previously, 
T) will always be negative and contribute to a repulsive force, no matter how much the 
anisotropy (M 3 ) is. And dF is the one that the anisotropy is mainly associated with. From 
Fig.3 we can see that T 2 takes the dominant place in most of the range. It can be either 
positive or negative depending on the anisotropy. Although t\ > e 3x > e 2 always produces 
a repulsive force, the total force will not necessarily to be repulsive, if the anisotropy makes 














FIG. 4. T, and T 2 vs. M 3 . Mi = 1.5 and M 2 = 1.1. The two vertical dash-dot lines 
divide the range of M 3 into 3 parts corresponding to e 3z < e 2 , e 2 < e 3z < €\, and e 3z > 
form left to right, respectively. The positive value corresponds to the attractive force, 
while the negative value corresponds to the repulsive force. 


e 3z to be out of this range, as shown in Fig.3. This is a result of the competition between dq 
and T 2 . If the anisotropy makes T 2 to be positive and have a greater amplitude than 'F 1 , 
the direction of the force will switch to attractive from repulsive. But without anisotropy, 
T 2 will be surely negative and the force will be always repulsive. To achieve an attractive 
force, anisotropy is a must in this case. Therefore, we can conclude that it is the anisotropy 
of the media between the plates that produces the attractive force in this case. 

This can also be understood in an intelligible and vivid manner as flowing. According to 
Lifshitz’s theory[3], the attractive force arises when e 3 is out of the range (e 2 ,ei). However, 
in above case, e 3x is assumed to be just in (e 2 , ei), which may always contribute to repulsive 
force. To produce an attractive force, one must let e 3z to be a little far form this range. So 
that the ” average value” (the quotation marks here mean that it is not the real mathematical 
average value 19]) of the e 3 can be out of this range, and the attractive force can arise. And 
to make e 3z out of the range (e 2 ,ei), anisotropy is a must in this situation. This is why the 
attractive force (T > 0) happens only on the left and right edges of Fig.3. 
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Now let’s turn to the case that e 3a: is not in the range (e 2 ,£i). If we let £32 < 62 (Mi > 
M 2 > 1), 'hi will be always positive and produce an attractive force. As shown in Fig. 
4, 'ho which is associated with anisotropy is still dominant in most of the range. Without 
anisotropy, \& 2 will be positive and the force will be always attractive. To achieve a repulsive 
force, M 3 must not be 0 in this case. We can also conclude that it is the anisotropy of 
the media between the plates that produces the repulsive force in this case. According to 
Lifshitz’s theory |s|, the force is repulsive only when e\ > e 3 > e 2 (or e 3 < e 3 < £ 2 ) is 
satisfied, and it is attractive in all the other cases. However, it can be found in Fig.4 that 
the total force can be repulsive, even when both e 3x and e 3z are out of the range (£ 2 , 61 ) ! 
The explanation is similar to the previous case. As shown in Fig.4, if e 33 . < e 2 , to achieve 
a repulsive force, M 3 (e 3 ~) must be larger than M 2 (£ 2 ) to bring the ’’average value” of e 3 
in (£ 2 ,£i). However, it should not be too much greater than M\ (ei), as this may make the 
’’average value” of e 3 beyond t\. That is why the repulsive force only appears in the middle 
region of Fig.4. 

From above discussion, we can find that the direction of the force is affected by M 3 which 
refers to the anisotropy of the media between the plates. If anisotropy of the in-between 
material makes the ’’average value” of e 3 to be in the range (e 2 , £ 1 ), the force will be repulsive. 
Otherwise, the force will be attractive. The border curve defined by T = 0 is shown in Fig.5 
(for convenient, M\ is fixed to be 1.5). The regions with T > 0 are the attractive regions, 
while the regions with T < 0 are the repulsive regions. 

IV. CONCLUSION 

We have calculated the Casimir force between two plates separated by a gap filled with 
anisotropic media. The result shows that the anisotropy of the media between the plates 
plays an important, or sometimes dominant, role in the total force. The more important 
thing is that it can affect the direction of the force. If ignoring the anisotropy of the in- 
between material makes the force to be repulsive, by contrast, taking the anisotropy into 
account can produce an extra attractive force; and if ignoring the anisotropy of the in- 
between material makes the force to be attractive, taking the anisotropy into account can 
produce an extra repulsive force. This can be explained as the result of the competition of the 
attractive term and repulsive term in the force function. And it can also be understood, in an 
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FIG. 5. The repulsive and attractive regions in M 2 — M 3 plane. The solid curve 
corresponds to T = 0 = 1.5) 

intelligible and vivid way, as the anisotropy makes the ’’average value” of the permittivity of 
the intermediate material to be in a certain range to produce the force in the corresponding 
direction. Finally, we should mention that it is still not easy to observe the above effects 
experimentally, since it is difficult to achieve high anisotropy in liquid. 
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Appendix A: The detailed derivation of the Casimir energy 


In this section we will present the detailed derivation of Eqs. (]6]) - (l8j) . The electric and 
magnetic held polarization vectors e q and h^ in Eqs. tHjl and ([5]) with frequency ou can be 
determined according to the classical Maxwell equations (In SI Units) 


d 

V x E = B 

at 

V x B = 44 e *E 
c z at 


(Al) 

(A2) 


€i (*=1,2,3) is the relative permittivity in region i, and its detailed representation can be 
found in Eqs. (H])-(J3]) , respectively. The solution of the electric and magnetic fields can be 
expressed in the form of plane wave as 


E 

B 



e* 

hr 

by 

\ b j 


(otx-\-'yz)—iujt 


(A3) 


e x , e y , and e z are the elements of the electric held (E) polarization vectors, while b x , b y , and 
b z are the elements of the magnetic held (B) polarization vectors. As we have assumed that 
all the materials considered in our work are nonmagnetic, we can have b x = Hoh x , b y = /J.oh y , 
and b z = /U 0 /r 2 . 

In Region I (as shown in Fig.2), substituting Eq. (JT|) and Eq. (IA3D into Eqs. (IA1D and 
(1A2D . we can have the eigenequations of the transverse elements of the electromagnetic held 
polarization vectors e x , e y , b x , and b y . 


/ 0 0 l-^o\ 

£1 


f e x ^ 


* e-x ^ 

0 0 0 -1 


Gy 


Gy 




= 7 


d 0 0 0 


Cby 


Cby 

y 0 —6i T o? 0 0 ) 


\ cb x / 


\ cb x J 


The eigenvectors and the corresponding eigenvalues can be written in the matrix form as 
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following: 


W 1 


( 0 0 



0 f 0 

1 1 

10 1 

0 10 / 


V = (7i 72 73 74 ) = ( -i^i -i^i ' lt i i*i ) 


(A5) 


(A6) 


with t\ = a 2 — ei- Each column of the matrix W 1 represents an eigenvector. Four eigenvalues 
correspond to four independent mode solutions. The general solution of the electromagnetic 
held should be the linear superposition of the four mode solutions, and the superposition 
coefficients are the amplitudes of each mode. Then the transverse elements of the electro¬ 
magnetic held in region I can be written as 


( E x ^ 


( /^K 0 (ax+~t\z) \ 

Ey 

= W 1 

j^^iKo{ax+^\z) 

CBy 


J^QiK 0 {ax+'i\z) 

\ cB x ) 


y A^QiK 0 (ax+'y\z) j 


(A7) 


where A* (j=1,2,3,4) is the amplitude of the jth mode solution. 

Similarly, in region II, where it is also isotropic, the transverse elements of the electro¬ 
magnetic held can be written as 


with 


and 


( E x ^ 


( J^IgKo(ax+~iY z ) ^ 

Ey 

= IE 11 

^II e iA'o(aa:+72 I z) 

CBy 


^II e i/Co(aa;+73 I 2 : ) 

\ cB x ) 




W u 


( 0 

_i_ 

t2 

0 

V 1 


j*2 Q i*2 \ 

£2 £2 

0 f 0 

£2 

1 0 1 

0 10 / 


7* = ( 7i : 7? 7? 74 ) = ( —i^2 -i^2 i^2 ^2 ) 


(AS) 


(A9) 


(A10) 


where t\ = a 2 — 62 , and A * 1 (j=l,2,3,4) is the amplitude of the jth mode solution. 
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In region III, where it is anisotropic, the case is different, as the permittivity has different 
form as shown in Eq.flHJ). Substitute Eq.(j3]) and Eq. (1A3[) into Eqs. flAlli and (IA2D . and we can 
have the eigenequations of the transverse elements of the electromagnetic held polarization 


and by. 





( 0 

0 

1 — — 

0 N 


e 3z 


0 

0 

0 

-1 


£ 3x 

0 

0 

0 


\ 0 

—e^x + « 2 

0 

0 ) 


f e x \ 

e y 

Cby 

\ cbx ) 


= 7 


( e x \ 

e y 

Cby 

\ cb x ) 


(All) 


The eigenvectors and the corresponding eigenvalues can be written in the matrix form as 
following: 


IE 111 = 


/ o 0 

£3x t3x 


(A 12 ) 


0 10 1 

\ i oio/ 

7 m = (7? 1 7? 1 7? 1 7i n ) = ( -its, ) (A13) 

with t 2 x = a 2 — €s x and t\ z = ( a 2 — es z )(e^x/^z) ■ Similarly, in region III, the transverse 
elements of the electromagnetic held can be written as 


( E x \ 


( y^llQiK 0 {ax+^Y lz ) ^ 

Ey 

= IE 111 

e \I<o(ax+~t™ z) 

CBy 


j^Ul e iK 0 (ax+'vl u z) 

\ cB x ) 


^ J^ l Q lK o( ax +l™ z ) j 


—\Ujt 


(A 14) 


where A / 1 (j=l,2,3,4) is the amplitude of the jth mode solution. 

As the surface modes should be exponentially decaying for z > 0 and z < a, we have 
A 3 = = A\ ] = A% = 0 BQ. As the transverse elements of the electromagnetic held are 

continuous at z = 0 , we have 


/ 0 -iA \ 


1 1 


0 


0 1 
V 1 0 / 


A 

Al 


= W 


hi 


( A™ \ 
A” 1 
Af 

\ Ai H ) 


(A15) 
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As the transverse elements of the electromagnetic held are continuous at z = a, we have, 

/ 0 & \ 

C2 
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t2 

0 1 

V 1 0 ) 




Kot2a 


A^Q~Kot2d 


= w 


III 


A^Q^ot^zO, 


y^IHg— Kot.3 x a 

\ A l l l e~ Kot ^ a ) 


(A16) 


Eqs. (lA15|) and (1A16I) include eight linear homogeneous equations relating the unknown 
parameters A\, A\, Ag 1 , Aj, A” 1 , A™, Ag 11 , and A^ 1 . It has the nontrivial solutions if 
the determinant of its coefficients is equal to zero, which leads to the equation for the 
determination of the proper frequency u ; . And the equation of determinant of coefficients 


20] equalling zero can be transformed into following form 

Y * [Gh(o;j_) * G 2 (o;||)] = 0 (A17) 

where Y is a function that is not always equal to zero. The subscripts _L and || indicate 
the mode with the polarization of the electric held perpendicular and parallel to the plane 
formed by k\\ = (k x , k y ) and z (As we have assumed k\\ to be parallel to the x direction, 
this plane is also x-z plane). Introducing u> = i£, the functions and G , 2 (<^||) can be 


expressed as following (Using the re 


t 2 z into the functions of p or ,s'i i2 21|) 


ation p 2 = 1 — a 2 /e 3x , we can transform t 2 , t 2 x and 


n , \_ncc\- i (si-p)(s 2 -p) 

Gi{w±) - G i(i£) — 1 - -— w——v-exp 

(si +p)(s 2 +p) 


2 pat; 


\/ c 3 x (|^^8) 


(t3xSl — Cl 1 ' )+f3 ‘ ) (e3a,S 2 — e 2 \/^ £ ^ _ 1 ' )+C3 ‘ ) 2a £\/ e 3 a > 2-1) + e 3 2 : e 3x 

G 2 (oo ||) = G 2 ( iO = 1- V ; - V , ^ (S*9) 


( e 3z s l + e l \J Ci ^ P e ^ +ei ~ ) (e 3x S2 + C2 \j 


e3z(p 2 -l)+e3z 

£3 z 


Eq. flAlSD is the same as Eq.(!71). Substituting Mi = ei/e 3x , M 2 = e 2 /e 3x , M 3 = e 3z /e 3x and 



E 


U). 


n, || 


(A20) 


The summations over n can be performed with the help of the argument theorem which has 
been applied in Refs. { 4 , 9, 17]. And then we have 

h 


E{a) = 


87 r 3 


/*27T /»oo 

kdk dO (]nGi + lnG 2 )d£ 
Jo Jo 


(A21) 


Substituting a 2 = e 3x (l — p 2 ) and k = into Eq. (1A21l) , we can arrive at Eq. 
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Appendix B: The approximate calculation of the integration in the Casimir force 


In this part we will introduce how we approximately calculate the integration over £ in 
the Casimir force in the limiting case. The Casimir force in Eq.flHl) can be rewritten as 
following 

F = Fi + F 2 


with 


Fi = 
F 2 = 


h 


27r 2 c 3 
h 

2i r 2 c 3 


p 2 dp 


pPdp 


A/ 2 £3 


-3x 




3/2 1 3 


“3a; 




1-G 1 

~Gi 

'1-G 2 

Go 


(Bl) 

(B2) 

(B3) 


The detailed expressions of G\ and G- 2 can be found in Eq. (JTJ) and Eq. (J51). 

For the limiting case that the separation a is larger than the characteristic absorption 
wavelength of the material, we can replace e 1 , e 2 , e 3x and e 3z by their values at £ = 0, i.e. 
the statistic dielectric constants j^|. We introduce Eq. (IB2D a new variable of integration 
X = 2pa^^/e 3 ^/c, and the integration should be taken over X and p. Eq. (jB2|) becomes 


F\ 


he 


dp 


327r 2 a 4 1 /ei^ J l p 2 J 0 

sio +ps 20 + p , v \ 
-exp(A) 


X 3 dX 


-l 


_si 0 p s 2 o p 

The integration over X can be taken by using approximate equation f c 


00 x n dx 


n! JO mexp(x) — l 


3 he 

f°°dp 

"sio ~ps 20 ~p' 

I67 r 2 a 4 ^/e^ 

h P 2 

_s 10 + ps 20 +p_ 


F\ 

similarly, we can have the approximate result of F 2 

F, 3RC 


(B4) 

iQ- 

(B5) 


16ff 2 a 4 v/ei7 


f°° pdp 

Sio — MiP s 2 o — M 2 P 

A i >:i 

sio + M X P s 2 o + M 2 P 


(B6) 


Then we can have the result in Eq. (Eli. And from Eqs. (fl4l) and (fl5|l we can hnd that 

3 he 


F\ = 
F 2 = 


167r 2 a 4 yTi^ 
3 he 
167r 2 a 4 v / ei^ 


Ti 


T, 


(B7) 

(B8) 
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